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Abstract 

Relying on an exact time evolution scheme, we identify a novel transient energy transfer phe¬ 
nomenon in an exactly-solvable quantum microscopic model consisting of a three-level system 
coupled to two non-Markovian zero-temperature bosonic baths through two separable quantum 
channels. The dynamics of this model can be solved exactly using the quantum-state-diffusion 
equation formalism, demonstrating finite intervals of unidirectional energy flow across the system, 
typically, from the non-Markovian environment towards the more Markovian bath. Furthermore, 
when introducing a spatial asymmetry into the system, an analogue of the rectification effect is 
realized. In the long time limit, the dynamics arrives at a stationary state and the effects recede. 
Understanding temporal characteristics of directional energy flow will aid in designing microscopic 
energy transfer devices. 
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The analysis of simple-prototype quantum energy transfer problems [1-3] assists in eluci¬ 
dating fundamental thermodynamic concepts in open quantum systems [4]. In the common 
construction, energy flow through a quantum system is generated by coupling it to two 
macroscopic objects, thermal reservoirs of different temperatures [5, 6]. Alternatively, direc¬ 
tional energy flow can be attained by supplying work into an asymmetric system [7]. Energy 
transfer problems are interesting for exploring the foundation of classical statistical mechan¬ 
ics, quantum dynamics, and the crossover between the classical and quantum worlds [8, 9]. 
For example, understanding the emergence of the Fourier’s law of heat conduction from the 
principles of open quantum systems is a long standing problem [10-14], Moreover, achieving 
control over energy flow is of an enormous importance in many areas of science and tech¬ 
nology, including energy management in functional nanoscale devices [15-18], realization of 
information processing and computation in open quantum systems [19], control over molec¬ 
ular reactivity and dynamics [20], and refrigeration in metal-superconductor junctions [21], 

In microscopic devices, introducing a spatial asymmetry within an anharmonic struc¬ 
ture can result in different magnitudes for the forward and backward currents, under the 
application of a reversed temperature bias [22], This diode-like behavior had recently at¬ 
tracted considerable theoretical and experimental attention, including the demonstration of 
phononic [23], electronic [24-27] and photonic [22, 28-33] rectihcations. 

Given these developments, it is highly desirable to identify minimal conditions under 
which a diode-like behavior can be obtained, controlled, and enhanced [3, 20], Studies 
of rectihcation and unidirectional energy flow in quantum devices [2, 3, 16] were typically 
performed under certain-standard approximations (semiclassical operation, neglecting coher¬ 
ences in the subsystem, assuming a unique steady state), adopting quantum master equation 
approaches. Few works had considered design principles on exactly-solvable quantum models 
[34], This problem is fundamentally important: Can we derive, from microscopic quantum 
theories, sufficient or necessary boundary conditions for realizing a certain nonlinear energy 
transport in atomic, molecular or a nanoscale system? 

In this work, we aim at achieving unidirectional energy flow in an open quantum system. 
While our results correspond to a temporal behavior, they expose ingredients for asymmetric 
dissipation thus potentially, asymmetric nonlinear transport. We employ reservoirs with 
different (non-Markovian) spectral properties rather than with different temperatures. A 
structured-non-Markovian environment is characterized by the correlation timescale of its 
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fluctuations [35], while in a Markovian bath the memory time is shorter than any other 
characteristic timescale of the system of interest. The bath memory function dictates the 
manner in which information and energy flow from the system to the attached macroscopic 
bath, and the back-action of the bath on the system. A finite memory time is crucial for 
achieving control over the state of an open quantum system [36]. Essentially, there is nearly 
no revival of the system’s hdelity when it is attached to a memoryless Markovian bath [37]. 
Therefore, if the system is coupled to two baths of different memory functions (correlation 
times, as dehned below), a unidirectional flow can emerge: energy is fed back from the non- 
Markovian bath to the system, and simultaneously, the system is releasing its energy to the 
more Markovian bath. 

We investigate the dynamics of our model by employing a nonperturbative master equa¬ 
tion, derived from the quantum-state-diffusion (QSD) equation [38] (see Method). We 
show that due to the assignment of distinct memory properties to the baths, a transient 
unidirectional energy flow develops in a prototype model for energy flow across a quantum 
open system. The model consists a three-level system [39], one of the simplest realizations 
of quantum engines [40], and two uncorrelated bosonic baths, bulk objects. Furthermore, 
we manifest that we can control the magnitude of the temporal energy flow by introducing 
a spatial asymmetry into the system, coupling it with different strengths to the contacts. 
This asymmetry grants an effect which can be categorized as a “transient diode effect”: 
the magnitude of the energy flow is different under forward and reversed operations, upon 
interchanging the channels connecting to the reservoirs. 

Results 

Unidirectional energy flow. To investigate the flow of energy in our system, purely 
induced by the distinct environmental memory functions, the two baths are assumed to be 
at zero-temperature and the three-level system is assumed to be of a degenerate A-type, 
where the energy splitting between the high level |3) and the two lower levels, |2) and |1), 
is set as uj. The setup is shown in Fig. 1, where the energy current across the system is 
determined by the two energy flows between the system and the two baths. We use the 
following form for the correlation functions of the two baths, aj{t,s) = with 

j = 1,2, where Tj is the coupling strength of the system to the jth bath. This form 
corresponds to a Lorenz spectrum. When —)■ 0, the jth bath is eminently non-Markovian 
with long memory time as desired. On average, the dissipation rate of the system to the jth 
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FIG. 1: Schematic representation of our model, a degenerate A-type three-level system coupled 
to two zero-temperature uncorrelated baths characterized by different memory parameters 7 and 
coupling constants F. 

bath is greatly suppressed with decreasing 7 ^- [38]. In contrast, when 7 ^ 00 , the jth bath 

is memoryless. In this case, energy flow from the system to the bath is fully irreversible. 
Therefore, l/ 7 j could be used to measure the environmental memory time. 

As our initial condition, we only excite state-3, pssit = 0) = 1. We demonstrate now the 
transient energy transfer behavior in the system when ai{t, s) 7 ^ 02 (t, s), before reaching the 
stationary long time solution P 33 = 0. In Fig. 2(a) we use 71 = 0.2a; and 72 = lOo; and show 
that initially energy is released into both baths simultaneously. However, around ut = 2.2 
the system begins to absorb energy from bath- 1 . As presented in Method, energy flow 
then becomes unidirectional, directed from the left reservoir (bath- 1 ) towards the right side 
(bath-2), i.e., TZlFi] < 0 and 7l[F2] > 0. We end our simulation when P 33 becomes extremely 
small and pu and P 22 reach the stationary states. The region of this unidirectional flow is 
embedded within a green-dashed frame. 

The dynamics can be made more involved if both reservoirs are highly non-Markovian. 
In Fig. 2(b) we use 71 = 0.2a; and 72 = l.Oo;, resulting in multiple-alternating regions 
of bi-directional and unidirectional energy flow. Particularly, we observe three intervals 
(distinguished by the green-dashed frames) of unidirectional flow of energy, of nearly the 
same duration yet shrinking amplitude, before full relaxation of the excited state is reached. 
As expected, the directional flow takes place from the reservoir with longer memory time, to 
the side with shorter memory time, since at zero-temperature a completely Markovian bath 
can only absorb energy. The comparison of Figs. 2(a) to 2(b) also reveals that, as expected. 
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FIG. 2 : (Color online) (a) and (b) Population dynamics of the three states in the A system, as 
well as our measure for energy flow, TZ[Fj{t)\ with Fi = F2, 71= 0 . 2 a; and different values for 72: 
(a) 72 = 10 a;, (b) 72=0;. The L —)• i? notation indicates that the unidirectional flow proceeds from 
the left bath to the right one within the marked frames 

. (c) Duration time of unidirectional energy flow (in its first occurrence) as a function of the 
inverse memory time 71 and 72. The dark blue region corresponds to cases with vanishing 

unidirectional flow (zero duration). 

in the latter case the total evolution time towards the stationary solution is longer than in 
the first case. 

A finite difference between the memory parameters, I 71 — 72 I 7 ^ 0, is a necessary yet 
insufficient condition for the emergence of unidirectional flow of energy. This is shown in 
Fig. 2(c), where we display the time duration of unidirectional flow in the first interval. 
Note that energy is flowing in opposite directions (i? —)■ L or L —)■ i?) in the regions below 
and above the diagonal in Fig 2(c). Recall that in Fig. 2(b) we show that there may be 
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more than one occnrrence of unidirectional transfer in the overall dynamics. We find that 
to observe the effect, it is necessary to employ a reservoir with a long memory time, for 
example, 71 / 0 ; < 0.5, and a second reservoir with a shorter memory time, I 71 — 72 ]/^ > 0.5. 
Diode-like phenomenon. So far, we have demonstrated that within a certain time interval 
energy may flow in a unidirectional manner only due to differences in the memory capabilities 
of the two reservoirs. In Figs. 2(a), 2(b) and 2(c), we used Fi = F 2 , and it is obvious that 
the direction and magnitude of the flow will be fully reversed upon the interchange of the 
quantum channels indicated by Li and L 2 , which is equivalent to an interchange in the values 
7 i and 72 . However, distinct memory times can not induce a diode-like phenomenon, an 
asymmetry in the magnitude of the unidirectional energy flow under opposite “polarities”, 
as we explain next. 

We now show that the dynamics may be furthermore controlled by including an asymme¬ 
try in the coupling strengths of the three-level system to the baths, Fi 7 ^ F 2 . The resulting 
behavior corresponds to the thermal diode effect, see Figs. 3(a) and 3(b). We consider the 
following two setups: (i) a “forward” configuration with 71 > 72 and Fi < F 2 ; (ii) a reversed 
geometry, in which we exchange the values of the memory times but keep the interaction 
energy as in (i), thus 71 < 72 and F^’s hold. We emphasize that the condition 71 7 ^ 72 
allows a unidirectional flow of energy, and the spatial asymmetry |Fi — F 2 I 7 ^ 0 , provides the 
transient diode effect, yielding different magnitudes for energy flow in setups (i) and (ii). 

We first examine geometry (i) in Fig. 3(a). We use 71 /a; = 5 > 72 / 0 ; = 0.2 and Fi = F 2 / 2 . 
Energy flows unidirectionally towards the more Markovian bath-1, as we found before, during 
the so called “diode-like behavior interval” 2.5 < ut < 5. In geometry (ii) we employ 
72 /a; = 5 > 7 i/a; = 0.2 while keeping Fi = F 2 / 2 . Energy now flows (during almost the same 
interval) towards bath-2. In geometry (ii) the ratio of the energy flow towards the more 
Markovian bath of 7 = 5a; is larger than in case (i), given the stronger coupling to this bath. 
It is interesting to note that when we modify F we largely affect the flow of energy into the 
more Markovian bath, as compared to changes in flow to the highly non-Markovian bath. 
We could explain this phenomenon by noting that F, the coupling strength of the system to 
the reservoirs, is the only parameter which determines the rate of energy flow to a Markovian 
bath. In contrast, the effectiveness of flow to a non-Markovian reservoir is predominated by 
the memory time of the bath, characterizing how effective it is in dissipating excess energy. 

The duration of the diode-like behavior is plotted in Fig. 3(b) as a function of 71 and 
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FIG. 3: (Color online) (a) Transient unidirectional energy flow and a diode-like effect in two 
configurations: (i) (blue lines) 71 = 5a; > 72 = 0 . 2 a; and Fi = F 2/2 = a ;/2 
, and (ii) (red lines) 71 = 0.2a; < 72 = 5a; and Fi = F 2/2 = a;/2. (b) Diagram of the diode-like 
time (the first period of unidirectional occurrence) in the parameter space of 71 and 72 . 

7 i — 72 . The calculation is performed on a configuration similar to geometry (i) with 71 
taken always larger than 72 , providing a unidirectional flow from bath-2 to bath-1. We hnd 
that the interval of the diode effect (considering the hrst interval) is highly sensitive to the 
memory time of the non-Markovian bath, while the difference 72 — 71 has a weaker overall 
effect. It is shown that a long diode-effect time is attainable when 72 /cu < 0.3; the overall 
duration of the diode effect is hxed once (71 — 72 )/^; > 4. 

Discussion 

We considered an exactly solvable model with two dissipation channels directed towards 
two reservoirs. We demonstrated that when adopting non-Markovian baths with different 
memory properties, the transient energy flow can become unidirectional, typically flowing 
from the highly non-Markovian to the more Markovian bath. This is the case as long as 
Iti — 72 1 is larger than a certain threshold. Our analysis departs from the regular thermo¬ 
dynamic setup in which energy flow is driven by a temperature gradient across the system, 
to consider zero-temperature situations with non-Markovian baths. Moreover, we showed 
that the magnitude of the energy flow can be controlled to achieve an effect reminiscent of 
the diode phenomenon, by coupling the system to two contacts with different strengths. In 
conclusion, a sufficiently large difference in I 71 — 72 I yields the effect of transient unidirec¬ 
tional energy flow; a sufficiently large difference in I 71 — 72 I, with the additional condition 
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of an asymmetrical coupling strengths between system and baths, results in a diode-like 
phenomenon. 

Our 3-level system could be realized in the triplet ground electronic state of a nitrogen 
vacancy (NV) center [41], when cu is regarded as the zero field splitting. Transitions from the 
lower degenerate levels to the upper state can be selectively addressed via optical fields [42]. 
Besides optical relaxation, one of the two channels could be realized through the vibrations 
of the diamond lattice and the atoms comprising the NV center point defect. In this case, the 
photonic field would serve as the strong non-Markovian bath while the phononic environment 
would act as the more Markovian one. 

Transient unidirectional flow of energy can be achieved in open quantum systems prepared 
in a non-stationary state, when coupled to two different structured environments. This be¬ 
havior, obtained in our work by a microscopic quantum model without further assumptions 
and approximations, can be exploited for estimating the relative memory capabilities and 
non-Markovianity of competing baths, for constructing nonlinear quantum devices for the 
transport of energy, and for controlling unidirectional energy transfer, and potentially reac¬ 
tivity, in molecules [43]. This is an exotic phenomenon induced by the boundary condition 
in quantum mechanics. Finally, the principles governing the dynamics of the present system 
at zero-temperature could be employed for exploring the dynamics of a finite-temperature, 
driven three-level system, to study the combined role of anharmonicity, non-Markovianity, 
driving, and asymmetry on energy transport phenomena [44, 45]. 

Method 

The system takes a A-type configuration, with one excited state |3) and two degenerate lower 
levels, |1) and |2). The excited state may decay to either of these lower levels, and these 
dissipation processes (referred as “channels”) are directed by two different baths [46, 47]. 
The baths are set at zero-temperature. The total Hamiltonian (h. = 1) is given by 

hftot = i^|3) (3|-|- 'y ^ UjkO^jf^ajk ^ ^yT^{gkajkL^j (1) 

j = l,2;fc j = i, 2 -,k 

Here, uj is the energy splitting between the upper level and the lower two states. The 
system-environment coupling operators are represented by Lj = |j)(3|, where j = 1,2. 
They open up energy transfer channels of the system into the jth bath, (ojA;) is the 
creation (annihilation) operator for independent mode-fc in the jth bath, ^yT^gk stands for 
the coupling constant between the system operator Lj and the /cth mode in the reservoirs. 


In our design unidirectional-transient flow is achieved by utilizing reservoirs with different 
two-time correlation functions aj{t,s), defined below Eq. (2). Furthermore, when Ei 7 ^ r 2 , 
a process analogous to thermal rectification can be realized. 

The wavefunction |tl'(t)) denoted the solution to the Schrodinger equation with the total 
Hamiltonian (1) in the interaction picture with respect to iJenv = 

define \'ipt{z*)) = ( 2 :i|(z 2 |d/(t)) where \zj) = Iif\zjk) stands for the tensor product of the 
Bargmann coherent states for the environment modes. The exact QSD equation [38], for 
the stochastic wave-function |' 04 (z*)), is given as 

dt\Mz*)) = I -*w|3)(3| + ^ [LjZ*^ - E,-(t)|3)(3|] I (2) 

I J=l ,2 J 

Here k correlated processes, describing the stochastic influence 

of the jth bath. z*j^ are individual Gaussian-distributed complex random variables, whose 
ensemble average is defined as M[zjtZ*^ = aj{t,s) = operator 

Fj(t)|3)(3| includes the effect of the jth bath on the system dynamics. The function Fj(t) = 
j^dsaj{t, s)fj{t,s) satisfies 

s) = [iu + Fi{t) + F 2 {t)]fj{t, s), (3) 

with fj{t,t) = 1 and j = 1,2. The corresponding exact master equation for the reduced 
density matrix p(t) = M[\'ijjt{z*)){\'ipt{z*)\] can be constructed via the Novikov theorem [48], 

dtp = Ha;|3)(3|,p] + ^ ^^Fj{t)[Ljp, L]] + h.c.j . (4) 

i=i ,2 

Equation (4) immediately yields the solution of time-dependent populations of the three 
levels, 

Pssit) = 

P22{t) = [ [F2 {s) + F2{s)]pss{s)ds + p22{0), 

Jo 

Piiit) = [ [Fi{s) + F^{s)]p 33 {s)ds + pii{0), (5) 

Jo 

where Fj{t) = J^dsFj{s). Since the system is prepared in its excited state, ^ 33 ( 0 ) = 1 and 
P 22 ( 0 ) = Pii(O) = 0, the time-dependent energy current, defined positive when flowing from 
the system towards bath -1 and bath- 2 , is given by 2u7l[Fi{t)]p33{t) and 2 a; 7 ^[F 2 (t)]p 33 (t), 
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respectively, where 7^[-] stands for real part. These expressions identify the current as the 
population relaxation rate times the energy difference for the transition. For example, state- 
3 decays to state-1 by giving up energy through channel-1 to bath-1. Population decays from 
level-3 to level-1 thus directly relates to the amount of energy flowing from the system to the 
attached bath. The transitory energy transfer can therefore be measured by the real part 
of dimensionless coefficients Fj{t), j = 1,2, without invoking confusion, before pss vanishes. 

Placing bath-1 (2) at the left (right) side of the system, we now identify different transport 
situations: (A) when 7l[Fi(t)] > 0 and 7l[F2(t)] > 0, the three-level system is releasing 
energy to both sides; (B) the system is releasing energy to the left while absorbing energy 
from the right side when 7l[Fi(t)] > 0 and 7l[F2{t)] < 0; (C) energy flows towards the right 
bath in the opposite scenario, 7l[Fi{t)] < 0 and 7l[F2{t)] > 0; (D) the system is absorbing 
energy from both reservoirs at the same time if TZ[Fi{t)] < 0 and TZ[F 2 {t)] < 0. Manifesting 
the development of scenarios (B) and (C), with a (transient) unidirectional energy flow as 
indicated by TZ[Fi {t)]TZ[F 2 {t)] < 0, is the objective of our work. If the memory functions have 
the same spectral form, ai(t, s)/Pi = a 2 (t, s)/T 2 , then Fi{t)/Ti = F 2 {t)/T 2 from Eq. (3). 
As a result, 7l[Fi(t)] and 7l[F2{t)] acquire the same sign at all times, thus a unidirectional 
energy flow across the system cannot be realized, even when Pi ^ P 2 . Evidently, to achieve 
our goal we should employ reservoirs with distinct memory properties. 
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